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In this paper we calculate the transfer coefficients for evaporation and condensation of mixtures. 
We use the continuous profiles of various thermodynamic quantities through the interface, obtained 
in our previous works using the square gradient model. Furthermore we introduce the Gibbs surface 
and obtain the excess entropy production for a surface. Following the traditional non-equilibrium 
thermodynamic approach we introduce the surface transfer coefficients which we are able to de- 
termine from the continuous solution. The knowledge of these coefficients is important for many 
industrial applications which involve transport through a surface, such as for instance distillation. 
In our approach the values of the local resistivities in the liquid and the vapor phases are chosen on 
the basis of experimental values. In the interfacial region there are small peaks in these resistivities. 
Three amplitudes control the magnitude of these peaks. Possible values of these amplitudes are 
found by matching the diagonal transfer coefficients to values predicted by kinetic theory. Using 
these amplitudes we find that the value of the cross resistivities is 1-2 orders of magnitude higher 
then the one from kinetic theory. The results of both kinetic theory and molecular dynamics sim- 
f-j-^ . ulations support the existence of small peaks in the local resistivities in the interfacial region. The 

Q ' square gradient approach gives an independent way to determine the transfer coefficients for sur- 

I faces. The results indicate that kinetic theory underestimates the interfacial transfer coefficients in 

^ real fluids. 
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I. INTRODUCTION 



^ . In two earlier articles (T], Q we developed the general approach for the square gradient description of an interface 
between two phases in non-equiUbrium 7i-component mixtures. Using that approach it is possible to determine the 
continuous profiles of all variables through the interface during, for instance, evaporation and condensation. In 
this paper we will use these results to obtain the transfer coefficients for heat and mass transfer through the liquid 
vapor interface. The values of these transfer coefficients, or even their order of magnitude, is extremely important 
for industrial processes which involve evaporation and/or condensation of mixtures. Among these processes is, for 
QQ _ instance, distillation, when one needs to separate components with different volatilities. As this involves evaporation 
. and/or condensation repeatedly many times, it is very important to know the exact effect of the surface. Some values 
T-H ■ of the interfacial transfer coefficients may favor transport of one component, while the others may favor adsorption 
, of a component at the surface. Of particular interest are the values of the cross coefficients, which contribute to 
■ reversible transport, and which are in most descriptions neglected jsj. 
^\ ', A number of different methods have been used to obtain the surface transfer coefficients for one-component systems: 
' experiments 4, 5, 6, 7], molecular dynamic simulations [1, H, [Hlj Ell , kinetic theory |1^, ,14, 15, 16]. In a paper 
^ ] coauthored by one of us [l3| the interfacial transfer coefficients obtained from the gradient theory for a one-component 
\^ ' system were calculated and compared to the data in the above references. Even for one-component systems the 
database of interfacial transfer coefficients is poor and these data are pretty scattered. The situation is even worse for 
mixtures. There are only few experiments available for several systems a-t a very restrictive range of conditions, 
i.e. for instance, at infinite dilution. No molecular dynamic simulations are available yet. The only source of the 
values of interfacial coefficients is kinetic theory [l^, [ly| . This theory is most appropriate for short ran ge p otentials 
and low density gases. There is evidence from molecular dynamic simulations for longer range potentials 12] that the 
coupling transfer resistivities for liquid-vapor interfaces of real fluids are substantially larger than those predicted by 
kinetic theory. 

In the first paper [l| we discussed the balance equations and the Gibbs relation for the square gradient model of 
mixtures. The Gibbs relation enabled us to derive the entropy production. It followed that if one uses as thermody- 
namic forces the gradient of the inverse temperature, V(l/r), minus the gradients of the chemical potential differences 
with the nth component divided by the temperature, — V[(/ij — /z„)/T], for j = 1, — 1, and minus the gradient 
of the velocity gradient divided by the temperature, — (Vv)/T, that the conjugate fluxes are the total heat flux, Jg, 
the diffusion fluxes of the first n — 1 components relative to the barycentric frame of reference, 3j, and the viscous 
pressure tensor. Linear laws relating these forces and fluxes could then be given. Together with the balance equations 
it is then possible to calculate the proflles of all the variables. In the second paper Q we gave details of the numerical 
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solution procedure for stationary states. We defined excess densities for an arbitrary dividing surface and verified 
that the non-equihbrium surface, as described by these excess densities following Gibbs, is in local equilibrium. Our 
work [H, 0] extended the analysis given by Johannessen et al. [isl. [Toj. 

Given the validity of local equilibrium for the description in terms of the Gibbs excess densities, it is possible 
to develop a description using non-equilibrium thermodynamics as explained in the monograph by Kjelstrup and 
Bedeaux [3|. This is much easier than the continuous description. As we will verify in this paper the expression for 
the excess entropy production of a surface has the general form 

a.^^J^X, (I.l) 

i 

In this expression Ji are the heat and mass fluxes through the surface and Xi are the jumps in the intensive variables 
across the interface. In non-equilibrium one uses the finite jumps of the temperature and chemical potentials across 
the surface, which lead to a non-zero entropy productions in the interfacial region. These jumps become the driving 
forces for the heat and mass transport through the interface. Following the traditional approach of non-equilibrium 
thermodynamics we then write the linear force-flux relations. These expressions use the interfacial resistivities or 
transfer coefficients which are the key interest of this paper. 

Having the continuous profiles of thermodynamic quantities obtained from the non-equilibrium gradient model we 
are able to calculate these resistivities independently. This gives a way to determine the coefficients and therefore a 
possible source for comparison with future experiments and simulations. 

In Sec. [n] we derive the expression for the local entropy production for stationary states in the continuous descrip- 
tion. In Sec. |III| we discuss the properties of the excess quantities in three-dimensional space. We consider how the 
stationary state condition simplifies the non-equilibrium expressions. In Sec. |IVj we obtain the expression for the 
excess entropy production in an interfacial region. In Sec. [V] we give the force-flux relations and discuss the interfa- 
cial resistivity coefficients. We consider different sets of coefficients which are associated with different variables: gas- 
and liquid- side coefficients, as well as mass and molar coefficients. Any coefficient of one set is determined by the 
coefficients of the other set and equilibrium properties of pure bulk components. These sets are therefore equivalent. 
In Sec. |VIj we proceed to a particular two-component mixture and specify further details. In Sec. |VII| we discuss 
the different methods to obtain the resistivities from a non-equilibrium continuous solution. We give the results of 
our analysis in Sec. |VIII| . We analyze extensively different aspects of the problem and find the values of parameters, 
which make the interfacial resistivities obtained from the continuous solution of the gradient model to match kinetic 
theory. We give a discussion and concluding remarks in Sec. |IX| . 



II. LOCAL ENTROPY PRODUCTION. 



A. Gibbs-Duhem equation. 



Consider a two phase n-component mixture. Let T be the temperature field in this region, ipi = fii — fi„ be the 
chemical potential differences and p be the scalar pressure, which in case of a planar interface coincides with the 
parallel pressure p|| . Furthermore, let u, s, v be the mass specific internal energy, entropy and volume respectively, 
p = 1/v he the overall mass density and = pi/phe the mass fraction of the i-th component. The explicit expressions 
for these quantities are given in 0]. It was found that the Gibbs equation for such a two phase system is given by: 

i—l 

where v is the barycentric velocity, d/dt is a substantial time derivative and we use the summation convention over 
double Greek indices. is the tension tensor, which is given by 

(r t) - n I V f ^^'('^^ | ^P^^' *^ + 'j^ d^tir^d^jirj^ 

"'^ ' dXa dxn ^ ' V dXa dxg dXa dxg J ^ dXa dxg 

' i—l ' ' i,j — l 

It is non-zero in particular in the interfacial region, where the gradient variables V p and V^i are significant. All 
thermodynamic densities (except the entropy) have gradient contributions. Foe explicit expressions we refer to 
These densities are related by 

n-1 
i=l 



Substituting Eq. pi.SP into Eq. pi.ip we obtain 
This is the Gibbs-Duhem equation for a two- phase multi-component mixture. Note that since ipi = fii — fin: we have 

which is the usual contribution to the Gibbs-Duhem equation associated with the chemical potentials. 
For a stationary state the derivative d/dt — and Eq. (|II.4p takes the following form 



V \ s— + + -v^^\ =0 (116) 



dxp dx0 ^--^ dx0 dxa 
where aap = p Sap + "faf3 is the thermodynamic pressure tensor 



B. Entropy balance. 

The entropy balance equation is 



p^^-V-Js+a, (II.7) 
dt 

with the entropy flux = Js.tot — psv and the entropy production as. These were found to be 

- J.-V^ - E J.-V^ - (11.8b) 

2 — 1 ^ 



where 3q and are the heat and diffusion fluxes 



3q =Je— pve— pV— 7r-V = Je— J„j (e +pw) — TT-V 

Ji = Pi (vi - v) = J^, - £,i J,„ 



(11.9) 



where tt = tTq^ is the viscous pressure tensor. The energy flux Jg and the mass fluxes Jj, = PiVj and J,„ — p'v are 
convenient quantities since in stationary states 

V-Je = 0, V-J^, =0, V-J™ = (11.10) 

Furthermore, it follows from Eq. pi.7p that in stationary state 

as = V-Js,tot = V-J, +pv-Vs = V-(J^,tot + psv) (11.11) 

Using pi.6p and the conservation laws under stationary conditions, it is possible to show that 

T ^ T T dxa T ^ ' 

i—l 

The expression for the entropy production, used in Eq. pi.l2p contains dependent fluxes J^^ and 3m and thus force-flux 
relations cannot be obtained from it directly. 
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III. EXCESSES IN THREE-DIMENSIONAL SPACE. 
A. Definition of an excess. 

The definition of an excess requires the normal direction n to be defined in the interfacial region. The surface 
may be curved and we may introduce curvilinear orthogonal coordinates (xi, X3) with rj_ = xi being the normal 
coordinate and r|| = (0:2, X3) being the tangential coordinates. 

Let x^'** and x^'* be the boundaries of the interfacial region at the gas and liquid side respectively. Let be a 
function defined in the surface region. Furthermore let cj)^, where superscript b stands either for £ or for g, be the 
function (j) extrapolated from the bulk to the surface region. The extrapolation is done using the description in 
homogeneous phases which does not contain gradient contributions. Outside of the interfacial region (j)^ and (/) are 
identical but inside the surface, (j)^ in general differs from (p. We note the following identity for the extrapolated 
functions 

/(x'''^r||) = 0(:r^■^r||) (IILl) 

Furthermore, for any function F 

(IIL2) 

in the interfacial region, since for the bulk functions outside the interfacial region it is identity. We note however, that 
even though Eq. (jllLip and Eq. pn.2p are exact, any numerical procedure will break these equalities. This happens 
because the extrapolation procedure usually involves polynomials in order to fit an actual curve, which introduces a 
non-zero error in the extrapolated curve. 

We then can define the excess (/((a:'*, r||) of a density 0(r) per unit of volume in the 3-dimcnsional space as^ [l^l 

1 '■^ 



</)(a:^r||) = — — / dxihih2h3r^{r;xn (111.3) 

"2 '''3 Jx^'" 

where 

</>^^(r; x") = (j){r) - ^^(r) e{x' - xi) - /(r) 6(2:1 - a;") (IIL4) 

Furthermore, hi = hi{xi,r^^) are Lame coefficients for curvilinear coordinates and hf = hi{x'',r^^). Given that 0(r) is 

a density per unit of volume, excess (l){x^ , r||) is a density per unit of surface. The excess depends on the position of 
the dividing surface x^, which is the coordinate of the surface in the normal direction, and the position r|| along the 
surface. 

B. Stationary state of a surface. 

Consider the entropy production given in Eq. (jIL12p . All the terms but last one have the form J-V0, where 
according to Eq. pL10[) V-J = and (p is some scalar function. Thus, J-Vi/) — V-(J0). We show in Appendix |A 2| 
that 

= {J^cbY - {J^d^Y + V^||0) (IIL5) 

where all the functions on the right hand side are evaluated at r**. 

For each flux in Eq. (jn.lO[l we can write Vj_ + V|| • J|| = 0. This gives an approximate relation for the order of 
magnitude 

^ . ^ (nL6) 

iS.xi_ Axil 



^ In the literature one also uses an alternative definition, which we show to be wrong in the Appendix lA II . 
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As was discussed in the interfacial region breaks the 3-dimensional isotropy of the system. In addition to a 
typical macroscopic size of the problem £, there exists the microscopic size S, the sm'face width, which is of the order 
of few nanometers. There are quantities which change drastically on the distances of the order S in the direction 
perpendicular to the surface. However, the significant change of any quantity along the surface may happen only on 
a length scale £, which is of the order of either radii of curvature or the system size. Because of this property of a 
surface, we may not expect the change of the parallel component of a flux on a macroscopic scale along the surface 
to be much larger then the change of the perpendicular component of that flux on a microscopic scale through the 
surface. For the fluxes for which changes |A^Jx| and |A|| J||| are of the same order of magnitude, Eq. pil.6p takes 
the form AJ/S w AJ/£, which can hold only if A J — 0, since S £. This means that both Aj_J± ~ and A|| J|| — 0. 
If |Aj^Jx| ^ l^ll"^!!! this statement becomes even stronger. We then may require that for a thin surface^ 

ViJi(r)=0 (III.7a) 
V||-J||(r)=0 (III.7b) 

Thus, a stationary state condition V- J = has a form of Eq. (|III.7p in an interfacial region. The extrapolated fluxes 
satisfy the same equation 

V±Ji(r)=0 (111.8a) 



V|i-J|;(r)=0 (111.8b) 



since the extrapolated flux fields also satisfy V-J*" = 0. 

Both Eq. (|III.7aP and Eq. (|III.8aP are first order ordinary differential equations which depend on a constant. These 
constants must be the same, since according to Eq. (jIII.l[) J^{r^''^) = J±{r^''') at boundary points. It means that 
and J± are the same functions^: 

Ji(r) = Ji(r) = Jx(r) (III.9) 
Consider the last term in Eq. piI.Sp . Since parallel divergences of both J|| and jj] are zero 

vi^cj)) = JirV|['/' (III. 10) 

where we used that V||6( xi) = 0. 

Substituting Eq. (|IlL9|l and Eq. (jllLlOp into Eq. pIL5|) we obtain 

= Ji(0^ - (t)<^) + J^V^^ (III.ll) 
where all the functions on the right hand side are evaluated at r**. 

IV. EXCESS ENTROPY PRODUCTION. 

Applying Eq. (jIII.ll[) to each term in Eq. (jII.12[) for the entropy production we obtain the general form of the 
excess entropy production for a surface in stationary state 

'^s — J^^tot^ ~ Js,totj_ + V|| ■ Js,tot|| = 



-^e^-L ( j ^'^i^^y-pl T9 J \ TS J (IV.l) 




( 



2 



For the special case of a system with planar surface in cartesian coordinated with all the fluxes directed perpendicular to the surface, 
which is the case studied in this paper, these equations follow straightforwardly. 

Note that Eq. I IIII.9II does not lead to the relation J^(x''''', ry) = Jj_{x^'" ,r^^). Eq. IIIII.9I I is the relation between values of different 
functions at the same point but not the relation between values of the same function at different points. However, it follows from 
Eq. IIIII.Tal l that in curvilinear coordinates d{h2h^J {dx^) = and therefore /i2/i3J± = const but not Jx = const. 
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where p,i = /li + v^/2 — g-r*. 

The next step of the analysis is to provide constitutive relations in order to relate thermodynamic forces 
to thermodynamic fluxes Jk for the whole surface. This requires that the excess entropy production has a form 
as =J2 JkXk- However, as one can see from Eq. (jIV.ip . the terms related to fluxes along the surface do not have this 
form. One has to make further assumptions on the nature of these terms to write them in this form. As our work 
is focused on transport into and through the surface we will not consider non-equilibrium perturbations which are 
applied along the surface. This guarantees that all the terms along the surface are equal to zero. The only nonzero 
component of any flux J is therefore Jj^, which we will denote simply as J. We will furthermore restrict ourself to 
non- viscous fluids. The expression for the excess entropy production simplifies to the following 

i\ " / ~t ~g 



£'4^-^ (.V-2) 



' 1=1 ^ 

It is convenient to write the excess entropy production in terms of the measurable heat flux J^, rather then the 
total energy flux Jg, which is defined as 

n n 
i=l 1=1 

where we used Eq. (|II.9p and hi = /i^ + v^/2 — g-r* = fli + Ts,;, where Si is the partial entropy and hi is the partial 
enthalpy. While Eq. (jIII.9|) is valid for Je and J^^ it is not valid for J^: the difference between the measurable heat 
fluxes extrapolated from the gas and the liquid side is 

n 
1=1 

In terms of measurable heat fluxes the expression for the entropy production becomes 



' f4 - ^) - E ^ (Af - + 4{T' T^)) (IV.5a) 
' i=i 



= ' f ^ - ^) - E ^ i^' - + 4iT' - T^)) (IV.5b) 



J'i / Zl^ J'g 

' i=l 

It is important to realize that Eq. (jlV.SP are exactly equivalent to Eq. (jIV.2|l . It is common to do these transformations 
neglecting third and higher order contributions in the deviation from equilibrium. Such approximations were not 
needed here. 

Eq. (jIV.SP has the form of the entropy production for the surface used in [3]. It was obtained there using the 
local equilibrium hypothesis, which we have proven to be valid in 2]. In this article we have derived Eq. (jlV.Sp 
independently, by calculating the excess of the continuous entropy production in the gradient model. 

V. SURFACE TRANSFER COEFFICIENTS. 

Consider Eq. (jIV.5[) for excess entropy production which has the form 

n 

= J,%-^Je.X, (V.l) 

i=l 

We will use the form (jV.ip further, specifying the explicit expressions for fluxes and forces where needed. Following 
the common procedure we write the linear force-flux relations for a given entropy production: 

(V.2) 

-Xj = RjqiTeq, Ipeq) Jq + YJl=l Rji{Teq , Ipeq) 
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As these relations are true only to the linear order in perturbations, we must keep only linear contributions in all 
terms. It means that all the resistivities in Eq. (|V.2p are functions of only equilibrium temperature Tgg and equilibrium 
chemical potential difference V'eg around which a particular perturbation is performed. They do not depend on the 
nature of the perturbation. 

Consider the following matrix notation of the above quantities 



X = 







^Rqq 


Rql . 


■ Rqn ^ 










Rlq 


Rll ■ 


■ Rln 








, R = 














Rnl ■ 


• Rn n J 







(V.3) 



Let fi indicate a measure of a non-equilibrium perturbation, so that X = X(6) and J = J(fi). Then Eq. (|V.2p can be 
written in a matrix form as 



X(6) =R(re„^e,)-J(fi) 



(V.4) 



For big values of fi Eq. (|V.4p is not correct, since big perturbations are not described by the linear theory. As we 
decrease fi, the accuracy of Eq. ()V.4[1 increases and in the limit fi ^ becomes exact. It means that Eq. (|V.4p should 
be understood as 



lim^X(fi) =R(re„Veg 



6^0 



•lim J(fi) 



(V.5) 



One should not write Eq. (IV. Sp in the form X(0) = R(req, -^eg) • J(0) however, as both X(0) and J(0) contain only 
zeroes and such an expression makes no sense. Even though X(fi) and J(fi) are continuous functions of 6, one should 
write limi3^oX(fi) and limg^o J(6) instead of X(0) and J(0) respectively. In practice there exists a particular measure 
fieg of a perturbation, such that for all fi < fieg Eq. (|V.4p is satisfied with a satisfactory accuracy. 



One should also note that the accuracy of a particular numerical procedure may limit the validity of Eq. (jV.5P as 
well. All the non-equilibrium profiles and therefore forces and fluxes are calculated by solving the system of differential 
equations numerically with some particular accuracy. If a perturbation rate is lower then this accuracy, say finum, 
the data obtained from the numerical procedure are not reliable. Performing a numerical analysis we must therefore 
replace the limiting value by fi„„m in Eq. (IV. 5p . 

As shown by Onsager [IH, the cross coefficients must be the same. We therefore have 



Rq 

R. 



Rii 
Ri 



(V.6) 



A. Gas- and liquid- side transport coefficients. 

For each of Eq. pV.5p one might associate the forces with 

Xq — -^q — ^ 

^ ^j^{i^-f^ + ^^(T'-T^)) (V.7) 

As we are in the context of the linear theory, however, we must linearize these forces with respect to the perturbation 
and discard all higher order terms. Leaving them would not increase accuracy but may affect the consistency of the 
linear theory. We therefore get the following phenomenological relations up to the linear order 

R'qqJq^' +EtlKJi' 

(V.8) 



1 1 _ 

eq 
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and 



(V.9) 



^ eg 



Here and after we omit arguments (Teq^ipeq) as long as it does not lead to confusion. 

The measurable heat fluxes are related by Eq. pV.4p which after linearization takes the following form 

n 

Jq'-Jq' = T.Ji^(!'U-hUq) (V-IO) 

1 = 1 

Comparing Eq. (jV.Sp and Eq. (jV.Op and using Eq. (|V.10[) we get the following relations between the coefficients 
associated with the gas an liquid measurable heat fluxes to linear order 







^qq 


^qq 






- 


j^e 

eq qq 


^^qi 


RS 

"'I, eq ''qq 




I- 


pe 

eq qq 


^R^ - 


RS 

"•'I, eq ^'-qq 


Rqi + K, 




pi 
eq qq 


^« 


R^ 

'\, eq ^jq 



(v-u) 

"'i.eq^^qq ^Hq '^i.eq^^qq 



, eq J, eq "qq 



where we took into account that /if_ = /if and /if, - /if = /if, e? " ^^1 eg- 

The coefficients on the one side determine uniquely the coefficients on the other side, having the values of jumps 
across the surface of the extrapolated enthalpies. It follows from Eq. (jV.lip that 



^qq 


''qq 








Rqi 


= R^ - 


iKeq 


-K 


) Rs 

eq) qq 


R-iq 


= R^ - 

^ Sg 


Keq 


-K 


eq) qq 




= R^. - 


Keq 


-hi 


eg) Rjq ~ ^hj, eg ~ h^ g^) i?f^ + (/if, — /if, g^) (/if, — g^) R^^ 



(V.12) 



One can notice, that symmetry of R^ coefficients implies the same symmetry of R coefficients and vice versa. 

To avoid confusion we recall that for the temperature T, chemical potential fij and the partial entropy Sj the 
superscript g or i means the value of the corresponding function extrapolated from either gas or liquid to the interfacial 
region and evaluated at a particular dividing surface x". We do not indicate which dividing surface is used as it is 
irrelevant for the present analysis. In contrast, J^'f = Jg — X^iLi hi{x^'^) J^^ where the partial enthalpy is evaluated 
the the gas-surface boundary x"'*. To the linear order however hi{x^''') ~ hi^eq{x^'^) — hf ^q{x'') = hf and since Jg 
and are constants, J^^ may be considered as the flux, evaluated at a dividing surface . Furthermore, superscript 
g or £ for the resistivity R^ or R^ neither indicate any position nor the extrapolated resistivity coefficient. It indicates 
the measurable heat flux with which the given resistivity coefficient is associated, either J^''^ or J^'^ respectively. 



B. Mass and molar transport coefficients. 



In applications it is common to use the mass flux of the components and the partial molar thermodynamic quantities, 
like, for instance, partial molar entropy. The above equations should use either molar fluxes and partial molar 
thermodynamic quantities or mass fluxes and partial mass thermodynamic quantities. The transport coefficients are 
different for different choices. Consider Eq. (jV.ip for excess entropy production. The thermodynamic forces Xi depend 
on partial thermodynamic quantities. We introduce therefore X™ as a force which uses partial mass quantities and 
as a force which uses partial molar quantity. Furthermore, let J^. = piVi and Jt^. = c^v^ be the mass and molar 
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flux respectively of the i-th component. As X™ = /Mi and J^. = J^.Mi, where Mi is the molar mass of the i-th 
component, the excess entropy production becomes 

n n 

a, = j;X, - J2 Ji. Xr = - ^ Jc. (V. 13) 



The force-flux relations become 



and 



- ^qq -Jq + 2^i=l ^qi -'i, 

Xq — Rqq Jq + ^gi '^Ci 

Vf r^u Jf 1 T^u J 



(V.14a) 



(V.14b) 



where the corresponding superscript for the resistivity coefficient indicates the association with the mass or molar 
quantities. These transport coefficients are related in the following way 



(V.15) 



^qq 


Dm 

— riqq 


Rqi 


= M,R^^ 


^iq 


^M,Rr, 




= MjM, R'Jl 



VI. TWO COMPONENT MIXTURE. 



We now restrict ourselves to the mixture considered in Q]. We consider the mixture of cyclohexane (1st component) 
and n-hexane (2nd component) in a box with gravity directed along axes x from left to right. The gas phase is therefore 
in the left part of the box and the liquid is in the right part of the box. The surface is planar. 

We have the following expression for the excess entropy production 



j/^Xq- J^.Xl- J^,X2 (VI. 1) 



and force-flux relations (jV.4p with 






( Rqq 


Rql 


R, 




Rlq 


Rll 


R 




\R2q 


R2I 


R: 




iq2\ 

X = \-Xi\ , R = I i?n i?i2 L J = I J^, j (VI.2) 

\—X2j \R2q R2I R22) X'^iiJ 

We will also use an alternative expression for the excess entropy production, which uses the total mass flux Jm = 
+ and the flux of one of the components, say J^^ = J^: 

(Ts = JqXq — J^X^ — JmXm (VI. 3) 

where X^ = Xi — X2 and X^ = X2- The resulting force- flux relations (|V.4p have the following terms 

(VI.4) 

where the coefficients from Eq. (|VI.2p are related to the coefficients from Eq. (|VI.4p as 

Rss 

(VI.5) 





xA 








, R = 




{-xt 







Rql 


- Rqi 


Rqm 


Rll 


— —Rm^ - 


R^m H~ Rmm 


Rlq 


= R^q 4 


Rmq 


R22 


— Rjnm 




Rq2 


— Rqm 




R12 


— Rmm ~ 


R^m 


R2q 






R2I 


Rmm 


Rm^ 
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Having the numerical solution for a particular non-equilibrium stationary state we know all the fluxes J and forces 
X used in Eq. (|V.4p : the constant fluxes are obtained directly from the non-equilibrium solution and the extrapolated 
bulk proflles are obtained using the procedure described in [2] . On the other hand we know only the local resistivities 
but not the resistivities R of the whole surface. 

Wc now consider the inverse problem: to determine the transport coefficients for the whole surface having the 
non-equilibrium solution. As one can see, Eq. (|V.4|1 has 9 unknown resistivities^ and only 3 equations. It is therefore 
not possible to determine all the transport coefficients uniquely having only one stationary state solution. In order to 
incorporate more equations we need to consider other non-equilibrium stationary solutions which are independent of 
the previous. As the transport coefficients depend only on equilibrium unperturbed state but not on non-equilibrium 
perturbations, considering different perturbations around the same equilibrium state we will get missing data. We 
must ensure however, that a given perturbation is small enough to be described by linear-order equations. This would 
require for instance Eq. (|V.6[) and Eq. (|V.11[) to be true. There are more constraints to be fulfilled which will be 
discussed in Sec. |VIII| . 

The non-equilibrium solution uses following profiles for local resistivities (see Q for details) 

rqq{x) = + [r^qq - r^^) qn{x) + agq(r^, + r^^) qi{x) 

•Tqiix) = + Ki -r^i)go(a;) +agi(r^i +r^i)gi(a;) (VI.6) 
rii{x) = rfi + (rfi -rfJgo(a;) -f- rf J qi (a;) 

where 90(2;) and qi{x) are modulatory curves for resistivity profiles which depend only on density profiles and their first 
derivatives. For each resistivity profile and are the equilibrium coexistence resistivities of the gas and liquid phase 
respectively. Coefficients a^^, otqi, an control the size of the peak in the resistivity profiles in the interfacial region. 
The non-equilibrium stationary state depends on the values of these coefficients. The surface resistivity coefficients 
R will therefore depend on these coefficients as on parameters, R = R(agq, a^i, an), which we will investigate. 

VII. METHODS TO OBTAIN RESISTIVITIES. 

We determine the transport coefficients from three different methods: from a "perturbation cell" method^, from an 
experimental-like procedure and from kinetic theory. 

A. Perturbation cell. 

Consider a stationary state which is perturbed from equilibrium by setting the temperature of the liquid^ T{x^) = 
(1 -I- f3T)Teq, the pressure of the gas p{x3) = (1 -|- (3p)peq and the mole fraction of the liquid (^{x^) = (1 -I- /3c)Cig 
independently. The resulting non-equilibrium state is therefore a function of parameters (3: 

X(/?T,/?p,/?c) =R(reg,^eg)-J(/?T,/3p,/?c) (VILl) 

where X, J and R are given by Eq. (jVI.2p . Consider the following set of 8 independent non-equilibrium perturbations: 



X( 


p, 


/3, 


P) 


= R{Teq,1peq) 


J( 


P, 


P, 


P) 


X( 


p, 


-P, 


P) 


= R{Teq,i>eq) 


J( 


P, 


-P, 


P) 


x( 


-/?, 


/?, 


P) 


= R(Te,,-0e9) 


J( 


-P, 


P, 


P) 


x( 


-/3, 


-P, 


P) 


= R(Te,,-0eq) 


J( 


-P, 


-/?, 


P) 


x( 


P, 


/?, 


-P) 


= R(Teg,-0eq) 


J( 


P, 


P, 


-P) 


x( 


P, 


-P, 


-P) 


= R(Te,,-0eq) 


J( 


P, 


-/?, 


-P) 


x( 


-/3, 


/?, 


-P) 


= R.(Teg,'0e9) 


J( 


-P, 


P, 


-P) 


x( 


-/3, 


-P, 


-P) 


= R(Teg,-0e9) 


J( 


-P, 


-P. 


-P) 



(VII.2) 



Solving the inverse problem we have to proof the Onsager reciprocal relations rather then impose them. 
^ This method was first used by Johannessen et. al. in [l7j for one-component system. Here we discuss the grounds for the legitimacy of 

this procedure and generalize it to mixtures. 
® One should not confuse T{x^) with T*. The former is the actual temperature at a; = , i.e. at the box boundary on the liquid side. The 

latter is the temperature extrapolated from the liquid phase to the interfacial region and taken at x = x", i.e. at the dividing surface. 
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Consider now the 3x8 matrices X and Z which contain 8 column vectors X and J respectively for each non-equilibrium 
perturbation specified above. For these perturbations X = X(/3) and 5 — are the functions only on one parameter 
(3. It follows from Eq. (|VII.2|1 that 

X{p)=RiT,g,ij,g)-Z{P) (VIL3) 

As it was discussed in Sec. [V], for practical purposes these limits should be calculated as the values of corresponding 
matrices at very small but finite value of f3. KiTi^gjip^q) depends therefore on (3 and we will keep it as an argument. 
From Eq. (jVII.3[) we obtain 

R(Te„ ^e,; /?) = ■Z'^iP)) ■ {m ■Z'^Wy' (VIL4) 

where superscript means the matrix transpose and means the inverted matrix. 

We note, that in order to obtain the resistivity matrix R uniquely, it is sufficient in principle to impose any 3 
non-equilibrium perturbations which have sufficiently small perturbation parameters /3t, Pp and P,^. As each of Pt, 
Pp and P(^ goes to zero the resistivity matrix will go to R(Teg, ipeq) as fast as all Pt, Pp and Pq go to zero. The method 
presented above makes the resistivity matrix converge to R(Teq, i/jeg) as fast as P"^ goes to zero, however. This is 
achieved by using 8 perturbations at the "corners" of a three-dimensional "perturbation cell", so changing P to — /3 
would not change the "perturbation cell" and the resulting R. 

Because of using 8 perturbations instead of 3, there are 5 superfluous perturbations which make the system of 
equations (|VII.3[) to be overdetermined. Contracting both sides of Eq. (jVII.3|) with Z^ we actually average all the 
perturbations which are spread around T^q and if^eq in the least square sense. As the components of Z matrix are 
linearly independent, this guaranteers the matrix to be invertible. Thus, the inverse matrix {Z-Z'^)~^ exists 

and Eq. (|VII.4p is mathematically legitimate. In the numerical procedure the expression on the right hand side of 
Eq. (|VII.4p is obtained using Matlab matrix division X/Z- 



B. Experiment-like procedure. 

In experiments it is convenient to measure the corresponding coefficients by keeping zero mass fluxes through the 
system. It is also convenient to work with the total mass flux Jm and the flux of one of the components , rather 
then with fluxes of each component separately^, J^^ and J^^. 

Consider a stationary state which is perturbed from equilibrium by setting the temperature of the liquid T{x^) = 
(1 + PT)Teq- The second perturbation condition is either = or C,^{x^) = Cfg and we introduce the perturbation 
parameter which is in the former case and 1 in the latter one. The third perturbation condition is either J„i = 
or p(x^) = Pf.q and the corresponding perturbation parameter P^ is or 1 respectively. The resulting non-equilibrium 
state is therefore a function of 3 parameters: 

X(/3t,/?4,/3™) = R{T,,q,ij,q)-:i{PT,P^,Pm) (VII.5) 
where X, J and R are given by Eq. (|VI.4[) . Consider the following set of 3 independent non-equilibrium perturbations: 

X(/3,0,0) =R(re,,^eg)-J(AO,0) 

X(/3,l,0) =R(re„7/;e,)-J(/3,l,0) (VII.6) 
X(/3,l,l) =R(re„V;eg)-J(/3,l,l) 

From the first of Eq. (|VII.6P we find 

Rqq {Teq,%i^eq) — Xq, QO / Jq,00 

R^qiTeq,^eq) - 00 /J^oo (VII.7a) 

Rmq{TeqT '4'eq) — A^m.OO / Jq, QQ 

where we use subscripts /3j/3„ instead the functional dependence P^, Pm) for simplicity of notation. From the second 
of Eq. (|VlL6|) we find 

Rqi {Teq^'^eq) = {Xq, 10 " Rqq {Teq,i>eq) Jq iq) / J^.W 

{T,q, Veg) = {X^, 10 - Riq {Teq, ^eq) Jq, lo) / J^. 10 (VII.7b) 
Rm^{Teq,i^eq) — {Xm, 10 " Rmq{Teq , ^^eq) Jq, lo) /J^.IO 



One of the reasons for this is that it is hard to make only J^^ = 0, keeping Jj^ finite. 
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The values Xio and Jio are found directly from the calculations and the values of Rqq{Teq,ipeq), Riq{Teq,tpeq) and 
Rmq{Teq,ipeq) ^re those which are found in Eq. (jVII.7a[) . given the perturbation rate (3 is small enough. From the 
third of Eq. (|VIL6|) we find 

Rqm {Teqjlpeq) — (-'^g, 11 ~ Rqq (Teq, Ipeq) Jq, 11 — Rq^ {Teq, Ipeq) J^, ll) /Jm, 11 

R^miTeqj^eq) (-''^5,11 ^ R^q {Teq , i^eq) J n - R^^ {Teq , ^eq) J^, ll) /Jm, 11 (VII.Tc) 

RmmiT^q , 'ip gq) — (^m, 11 Rmq{Teq , Ipeq) Jq^ 11 Rm^(Teq , Ipeq) J^, ll) / Jm, 11 

Again, all the quantities on the right hand side of Eq. (|VIL7c[) are known and we therefore can find the remaining 
resistivities. 



C. Comparison to kinetic theory. 



According to p. 180] kinetic theory gives the following expressions for the surface transport coefficients for two 
component mixture 



where 



257r V <^i 



16 Wi 



Rfi'iT,^) = (i?r)2 i?o (T, V) 1 + 32 



^2 



1 / 1 



1 3 

TT 4 



Ro{T, V) EE 2-9/2 (i?T)-5/2(cf / + 4/./M^y 
w,{T,^/j) =A,/(Ai+A2) 



(VII.^ 



(VII.9) 



where R is universal gas constant, and 4 eq ^'''^ thermal conductivity and the gas coexistence concentration 
of i-th component respectively, cji is the condensation coefficient, which is parameter in this theory and Sij is the 
Kroneker symbol. 



VIII. RESULTS. 



Using the procedures described above we obtain different sets of transport coefficients R(r, "0), each of them as a 
function of equilibrium temperature and chemical potential difference. Let us use subscript pc for the resistivity matrix 
obtained from the "perturbation cell" method and ex for the resistivity matrix obtained from the "experiment-like" 
method. In each method we calculate the resistivities associated with the gas- and liquid- side measurable heat fluxes 
using Eq. (fO)) and Eq. (fO]) . 

Furthermore we will use subscript kin for the resistivity matrix obtained from kinetic theory, for which only the 
gas- side resistivities are available. We calculate the transport coefficients associated only with mass properties. The 
corresponding molar coefficients may be calculated using Eq. (jV.lSp . As a result we obtain the following sets of 
resistivities: R^^, Kp^, Rf^, Rg^., all of which depend on temperature and chemical potential difference as well as on 
parameters aqq, aiq, an. In addition we obtain R^j„ which depend on temperature and chemical potential difference 
as well as on condensation coefficients tri and tT2. We have the following constraints, which they must obey for each 
T and tp'- 

- i) the second law consistency; 

- ii) the cross coefficients of each R matrix must satisfy Onsager relations (|V.6|1 : 

- iii) the corresponding components of R^^ and Rp^ as well as Rf^ and R^^ must satisfy Eq. (jV.ll[) ; 

- iv) the corresponding components of R^,,, Rf^. and R|-„ as well as R^^, and Rf^ obtained at the same T and ip 
must be equal. 

We study the dependence of the different resistivity coefficients on aqq, aiq and an and on T and ip and their 
convergence for small /3. We determine the values of the parameters for which the above constraints are fulfilled. 
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A. Onsager reciprocal relations. 

In this subsection we investigate the values of parameters aqq, aiq, an for which the Onsager relations are fulfilled. 
This is done for a particular values of equilibrium temperature and chemical potential difference T^q — 330 K and tpg^q — 
700 J/mol. In Tables [HUT] we give the relative error in percent for the gas-side cross coefficients \{Rij ~ Rj^) / R^j\-100% 
as a function of /? for aqq — 0, aiq = 0, an = obtained by different methods. 

TABLE I: Relative error in percent for gas-side cross-coefficients ob- 
tained by "perturbation cell" method at Teq — 330 and tpeq ~ 700 for 
different /3 and for aqq = 0, Uiq — 0, an — 0. 



/3 


Rql 


Rq2 


Rl2 


2.0e-002 


8.963066 


35.863259 


34.908631 


2.0e-003 


0.273286 


0.369082 


19.683274 


2.0e-004 


0.011726 


0.007231 


1.909391 


2.0e-005 


0.066375 


0.071266 


2.336652 


2.0e-006 


4.963895 


8.128243 


5.843913 



TABLE II: Relative error in percent for gas-side cross-coefficients ob- 
tained by "experiment like" method at T^q = 330 and tpeg ~ 700 for 
different /3 and for Uqq — 0, aiq — 0, an — 0. 



/3 


Rql 


Rq2 


Rl2 


2.0e-002 


1.275105 


0.828600 


754.982200 


2.0e-003 


0.038759 


0.363715 


38.708981 


2.0e-004 


0.131868 


0.238584 


6.247648 


2.0e-005 


1.301483 


2.056102 


20.984734 


2.0e-006 


13.282959 


20.788752 


632.124504 



As one can see, /3 — 0.02 is really an extreme perturbation and the difference is rather large. When we decrease (3 to 
2e-4 the differences become small. As we further decrease /3 to 2e-6 the inaccuracy of the numerical solution become 
comparable to the size of the perturbation. We conclude that the values for /3 to 2e-4 are closest to the converged 
values and use them as such. 

In Tables |IIHIIV| we give the same data for the higher continuous resistivities with rather substantial peak, when 
aqq = 10, aiq = 10 and an — 10. As one can see, the Onsager relations are fulfilled there again best for /3 = 2e-4 

TABLE III: Relative error in percent for gas-side cross-coefBcients ob- 
tained by "perturbation cell" method at Teq — 330 and ^peq = 700 for 
different /3 and for aqq — 10, aiq — 10, an = 10. 



/3 


Rql 


Rq2 


Rl2 


2.0e-002 


71.515410 


78.166809 


23.572836 


2.0e-003 


0.745604 


0.896547 


0.317348 


2.0e-004 


0.012358 


0.012650 


0.001919 


2.0e-005 


0.012078 


0.007485 


0.005290 


2.0e-006 


0.713969 


1.124994 


0.022121 



TABLE IV: Relative error in percent for gas-side cross-coefficients ob- 
tained by "experiment like" method at Teq — 330 and tpeq = 700 for 
different /3 and for aqq = 10, aiq — 10, an = 10. 
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Rql 


Rq2 


Rl2 


2.0e-002 


4.225362 


2.559393 


12.259260 


2.0e-003 


0.443944 


0.256804 


1.091842 


2.0e-004 


0.068621 


0.019788 


0.093041 


2.0e-005 


0.269764 


0.407090 


0.008844 


2.0e-006 


2.717575 


4.149484 


2.025054 



The similar picture is observed for the liquid-side resistivities and we do not give those data here. 
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We may notice that the behavior of the resistivities with respect to /3 is independent on the behavior of the 
resistivities with respect to aqq, aiq and an. This is natural, as these parameters control the different aspects of the 
system: j3 controls the perturbation rate, while a's are adjustable parameters, which control the size of the peak in 
the continuous resistivities. 



B. Second law consistency. 

In this subsection we investigate the values of parameters a^q, aig, an for which the second law of thermodynamics 
are fulfilled. That is that the excess entropy production is positive and therefore the matrix of the resistivity coefficients 
is positive definite. This requires that the diagonal coefficients are positive and for each pair ql, q2 and 12 of the 
cross coefficients the expression 

DR,k = RuRkk - ^{R^k + i?fc,)' > (VIII.l) 

must be positive. 

In Table [Vj we give the diagonal coefficients and expression (IVIII.1|) for each pair of the cross coefficients as a 
hmction of Ugg for aiq = 0, an = and (3 — 2e-4 obtained by the "perturbation cell" method. In Tables |VmVII| we 
give the same quantities for other choices of a. 

TABLE V: 2nd law consistency for gas-side coefficients. The diagonal 
coefficients and the quantities defined by HVIII.1|I . Data are obtained by 
"perturbation cell" method at T^q = 330 and tp^q = 700 for difi'erent aqq 
and for (3 = 0.0002, aiq = 0, an = 0. 



aqq 


Rqq 


Rii 


R22 


DRql 


DRq2 


DR12 





7.05644e-015 


0.0754717 


-0.0919278 


2.13025e-015 


-2.59473e-015 


-0.0277518 


1 


3.36047e-012 


0.0937784 


-0.0741586 


1.26056e-012 


-9.9683e-013 


-0.0278179 


10 


3.35408e-011 


0.259425 


0.0851534 


3.48053e-011 


1.14244e-011 


0.0874467 



TABLE VI: 2nd law consistency for gas-side coeflicients. The diagonal 
coefficients and the quantities defined by (|VIII.1[I . Data are obtained by 
"perturbation cell" method at T^q = 330 and tp^q = 700 for different aiq 
and for f3 = 0.0002, aqq = 0, an = 0. 



aiq 


Rqq 


Rii 


R22 


DRql 


DRq2 


DR12 





7.05644e-015 


0.0754717 


-0.0919278 


2.13025e-015 


-2.59473e-015 


-0.0277518 


1 


7.05608e-015 


0.0746391 


-0.0910331 


2.10664e-015 


-2.56935e-015 


-0.0271785 


10 


7.05304e-015 


0.0670813 


-0.0828915 


1.89251e-015 


-2.33855e-015 


-0.0222419 



TABLE VII: 2nd law consistency for gas-side coeflicients. The diagonal 
coefficients and the quantities defined by (|VIII.l|l . Data are obtained by 
"perturbation cell" method at Teg = 330 and V'e<j = 700 for difi'erent an 
and for (3 = 0.0002, aqq = 0, aiq = 0. 



an 


Rqq 


Rii 


R22 


DRql 


DRq2 


DR12 





7.05644e-015 


0.0754717 


-0.0919278 


2.13025e-015 


-2.59473e-015 


-0.0277518 


1 


7.05717e-015 


0.370078 


0.265626 


1.04468e-014 


7.49827e-015 


0.381226 


10 


7.10378e-015 


3.02063 


3.48284 


8.58316e-014 


9.89654e-014 


-69.2846 



We see, that the required quantities become positive for rather big values of aqq. They almost do not depend on 
the value of aiq and they are positive for moderate values of parameter an. It is clear that finite values of aqq and 
an are needed to have a positive excess entropy production. 

All the above quantities almost do not depend on the value of /9 in the range [le-5, le-3]. The "experimental-like" 
procedure leads to almost the same values of all the quantities. The liquid-side coefficients reveal a similar behavior. 



15 



C. Gas- and liquid- coefficients. 



In this subsection we investigate the validity of Eq. (jV.lip . In Table [Vlllj we give the relative error in percent 
between the left hand side and the right hand side of Eq. (jV.lip . 



TABLE VIII: Relative error in percent for invariant expressions in 
Eq. (|V.11[) obtained by "perturbation cell" method at Teq = 330 and 
ipeq = 700 for P = 0.0002 and Q,q = 1, ai, = 1, an = 1. 





11 


22 


9l 


ig 


g2 


2g 


12 


21 


0.000000 


0.000002 


0.000085 


0.000001 


0.000389 


0.000001 


0.000389 


0.000060 


0.000003 



For instance, the ql quantity is equal to - h[ ,^q i?^^) - {R^^ - h^^^ Rlq)\/\Rqi - M,e9 i?^,|-100%. The other 
quantities are defined in the same way. These errors almost do not depend neither on the value of (3 in the range 
[le-5, le-3] nor on the values of Uqq, aiq, an. The "experimental- like" procedure leads to almost the same results. 



D. Comparison to Icinetic theory. 



In this subsection we investigate the values of parameters aqq, aig, an which makes the coefficients agree with the 
kinetic theory coefficients. We do it for (3 = 2e-4 as this perturbation rate gives the most accurate results. We again do 
this for temperature Teq — 330 K and chemical potential difference tpeq = 700 J/mol. The values of parameters, used 
for kinetic theory are the same, as we use in our calculations. Particularly, the heat conductivities are Ai =0.0140 
W/(m K) and A2 =0.0157 W/(m K), Mi = 84.162 g/mol and M2 = 86.178 g/mol. We compare here only the 
"perturbation cell" method with kinetic theory. 

We found that the variation of aiq from to 10 makes the diagonal coefficients vary about 1% and the cross 
coefficients vary not more then 5%. As the variation of aiq is quite substantial, the variation in the coefficients which 
it induces is negligible. We therefore take aiq = in all further analysis. 



For the above parameters i?, 



gg, kin 



= 2.96792 X 10 . We found that Rqq,pc is practically independent on an while 



it depends linearly on a^g, see Fig. [T]. One can see from the plot, that R. 



qq, kin 



= i?, 



qq, pc 



for 



»99 



9. 




FIG. 1: Dependence of Rqq on aqq obtained by "perturbation cell" method at Teq = 330 and ^peq = 700 for aiq = and 
Qii = 1. Rqq. kin IS drawn as a constant line. 



The diagonal coefficients Rn.pc and i?22,pc depend both on aqq and an. Since we have found the corresponding 
to kinetic theory value of aqq by mapping the Rqq coefficient, we will further investigate the dependence of Rn.pc 
and i?22,pc using this value of aqq and varying only an. The diagonal coefficients Rn,kin and i?22,fcin depend, in 
their turn, on the condensation coefficients ui and CT2 respectively. We plot this dependence in the same plot with 
the dependency of Ra.pc on an, see Fig. [5]. The dependence of Ru.pc on an is given by the dotted line with the 
values of an drawn on the bottom a;-axes. The dependence of Ru, kin on ai is given by the solid line with the values 
tJi drawn on the top x-ax.es. 

Consider a particular value Ru.q of the diagonal coefficient Ra, where i is either 1 or 2, which is indicated by a 



horizontal dashed line on a figure. To find the value of an for which Ri- 



Rii, we draw a perpendicular from the 
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FIG. 2: Dependence of Rn^pc and i?22,pc on an (dots, bottom axes) and Rn^kin and R22,kin on ai and a2 (curve, top axes), 
respectively. Data are obtained at T^q = 330 and tp^q = 700 for aqq — 9 and aiq = 0. 

point where it crosses the dotted hne to the bottom axes. To find the value of ai for which Ru^ kin — Ra. o we draw 
a perpendicular from the point where the horizontal dashed line crosses the solid line to the top axes. For instance, 
the value R22.0 — 1-1 corresponds to an — 3 and 1J2 = 0.62. The value an = 3, in its turn, gives i?ii.o = 1-1 which 
corresponds to ui = 0.54. 

One may start by specifying an, rather then Ru^q, to find ai and 172. Then we draw a perpendicular from the 
bottom axes until it crosses the dotted line, which gives the value Ru^o of Ru^pc- Given the value of Raskin to be 
the same, we find the value of Ui as described above. For the above example an = 3 corresponds to ai = 0.54 and 
CT2 = 0.62. We see, that we may not specify both <ti and 172 independently: they must have the values which both 
correspond to the same an. For similar components, like those we are interested in, ai and CT2 should not differ much 
from each other, and therefore an, a coefficient which is related to the diffusion of one component through the other, 
should reflect this difference. 

Having the diagonal coefficient mapped we have the parameters aqq and an defined uniquely (and taking into 
account that aiq has negligible effect), as well as CTi and (T2 for kinetic theory. We now compare the values of the 
cross coefficients given by "perturbation cell" method and kinetic theory. 



TABLE IX: Gas-side transport coefficients obtained from kinetic theory 
and by "perturbation cell" method at Teq = 330 and ^})eq = 700 for 
P = 0.0002. 



parameters 


Rqq 


Rii 


R22 


Rql 


Rq2 


R12 


(Ti = 0.54 
(72 = 0.62 


2.96792e-011 


1.11091 


1.09136 


3.82826e-007 


4.41483e-007 


0.0130511 


OLqq — 9 














Ctlq = 


3.01874e-011 


1.12461 


1.13991 


2.31477e-006 


2.27003e-006 


-0.816559 


an = 3 















One can see from Table |IXj that while the diagonal coefficients are the same^, the cross coefficients we find are 
between 1-2 orders of magnitude larger than those found by kinetic theory. R12 even has a different sign. 



One should not expect exact compatibility between kinetic theory, which is most appropriate for gases with short range potentials, and 
the gradient theory, which is most appropriate for fluids with long range potentials. The purpose of this comparison in not to determine 
the exact values of adjustable parameters, but to show that it is possible to match coefficients in the two theories and to show the 
typical values of the parameters. 
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E. Temperature and chemical potential difference dependence. 

In this subsection we investigate the dependence of the resistivity coefficients on the temperature and the chemical 
potential difference. On Fig. we plot the these dependencies for Rgg, Rqi and Rn coefficients obtained from 
kinetic theory and "perturbation cell" method for the range of temperatures [325, . . . , 335] and for the range of 
chemical potential differences [400, . . . , 1000]. 




T/K 

FIG. 3; Dependence of Rqq on T and tp obtained from kinetic theory for cti = 0.54 and a2 = 0.62 (plane) and by "perturbation 
cell" method for aqq = 9, aiq = and an = 3 (points). 




J / K 325 400 

\\f I J/mol 

FIG. 4: Dependence of Rq\ on T and ^j) obtained from kinetic theory for ai = 0.54 and (J2 = 0.62 (plane) and by "perturbation 
cell" method for aqq = 9, aiq = and an = 3 (points). 

The domain of T and V' is not big, so the dependence on them is linear, as expected. 
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1.3^ 




T/K 325 400 y/J/mol 



FIG. 5: Dependence of Rii on T and ip obtained from kinetic theory for ai = 0.54 and (T2 = 0.62 (plane) and by "perturbation 
ceU" method for Oqq = 9, aiq = and On = 3 (points). 



IX. DISCUSSION AND CONCLUSIONS. 



In this paper we have studied stationary transport of heat and mass through the hquid-vapor interface in a mixture, 
using the square gradient theory [l], Q • We derived an expression for the excess entropy production of a surface from 
the continuous description, which is identical to the one derived directly for the discrete description using the property 
of local equilibrium Q. This makes it possible to give the linear force-flux relations for this case. These relations 
involve the interfacial resistivities or transfer coefficients, which were the main focus of interest in this paper. Given 
the numerical solutions of the non-equilibrium gradient model we were able to calculate these coefficients directly for 
a two-component mixture. This gives an independent way to determine the interfacial resistivities. 

The main input parameters of the model, aside from the parameters in the equilibrium square gradient model, 
arc the local resistivity profiles used to calculate the continuous solution. There is not much theoretical information 
about the numerical value of these resistivities. In the vapor phase one can use kinetic theory. In the liquid phase it 
is most appropriate to use experimental values. There is no experimental information about the local resistivities in 
the interfacial region. As the local resistivities change in the surface from one bulk value to the other, it is natural 
to assume that they contain a contribution similar to the profile of the order parameter. There is also evidence from 
molecular dynamics simulations for one-component systems p^ l that there is a peak in the local resistivities in the 
surface. As we are in the framework of the gradient theory, it is naturally to assume that this peak is caused by a 
square gradient term, which is similar to the gradient contribution, which the Hclmholtz energy density has in the 
interfacial region, namely |Vpp. The amplitude of this peak is not given by any theory and is used as a parameter. 
We therefore get that each of three local resistivities for a two-component mixture has the form given in Eq. (jVI.6p . 
Thus we get three adjustable amplitudes, aqq, aiq and an, two of which are found to contribute significantly to the 
value of the transfer coefficients. 

In order to determine the typical values of the a's we need to compare our results with independently obtained 
resistivities. Unfortunately, not much experimental data are available for multi-component resistivities and, to the 
best of our knowledge, no data are available for our system. Furthermore, no molecular dynamic simulations are 
available for mixtures. The only reasonable source of comparison is kinetic theory, which gives the expressions for 
the interfacial resistivities or transfer coefficients [l^ [l^ . We therefore compare our results to kinetic theory. Having 
three adjustable parameters in the gradient theory, a^q, aiq and an, and two adjustable parameters in kinetic theory, 
the condensation parameters cti and a2, we are able to match three diagonal coefficients Rqq, Ru and R22- We found 
that Rqq does not really depend on aiq and an. This makes it possible to fit aqq using Rqq alone. For the values 
of the temperature and chemical potentials considered this gave Uqq ~ 9. We furthermore found that the interfacial 
resistivities did not really depend on aiq. We therefore took this amplitude equal to zero. In kinetic theory Ru 
and i?22 depend on the condensation coefhcients ci and cr2, respectively. Choosing an = 3 gives values for the 
condensation coefficients of 0.54 and 0.62. As the components considered are very similar it is to be expected that 
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these coefScients are close to each other. The values of a's obtained from the matching are such that the excess 
entropy production of the surface is positive, the second law is obeyed and the Onsager relations are valid. Having 
found the values of the a's from the diagonal transfer coefficients the value of the cross coefficients follows uniquely. 

We found that the values of the cross coefficients, obtained by our method are between 1-2 orders of magnitude larger 
than those found from kinetic theory. This confirms results from molecular dynamics simulations |.12j , where it was 
found that increasing the range of the attractive potential increased in particular the cross coefficients substantially 
above the values predicted by kinetic theory. This is an interesting result, indicating that kinetic theory underestimates 
the transfer coefficients for real fluids. This also indicates, that the effect of coupling will be important in the interfacial 
region. Experiments also confirm the importance of the cross coefficients [1, [7| . 

The effect of cross coefficients can be related to the measurable quantities, such as measurable heat of transfer 
q* = —Rqi/Rqq. This quantity can be associated both with gas and liquid phases in accordance to the corresponding 
heat fluxes. The difference g*'^ — q*'^ = —{Rqi — Rqi)/Rqq — — — ^leg) is equal to the difference of partial 
enthalpies between gas and liquid in equilibrium. This quantity is substantial, which indicates that q*'^ — q*'^ is 
also substantial. This, in turn, makes the difference between the cross coefficients on gas and liquid side to be non- 
vanishing. This gives a theoretical ground for the importance of coupling in the interfacial region. Experiments 0, 0] 
confirm the size and importance of the heat of transfer on the gas side. 

We did the comparison for one value of the temperature and chemical potential only. If one extends the analysis 
to a larger domain, one finds that the a's depend on the tem per ature and the chemical potential difference; we refer 
to [l3| in this context. The results of kinetic theory [l^, [l^ llSl. [iGj and molecular dynamics [11] both support the 
existence of a peak in the diagonal local resistivities and therefore the use of finite values for aqq and an. 
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APPENDIX A: EXCESS IN CURVILINEAR COORDINATES. 
1. On the definition of an excess quantity using curvilinear coordinates. 

One may think of an alternative definition of an excess quantity 

^a{x',X2,X3) = da;i/ii(/)^^(r;x") (A.l) 

We note however, that 0q has no physical meaning, while (j) has. The reason for this is that 

<I> = // dx2dx3h2h^(t)ix'',X2,X3) ^ // / dxidx2dx3hih2h3(j)'''^{r;x'') (A. 2) 



S JJSJX- 



is the total amount of some physical quantity in the volume which is limited by the surfaces S at x^'^ and x^''', while 
$a = 11 dx2dx3hlhl'(f)a{x'' ,X2,X3) = ll / dxi dx2 dx^ hi (j)"'' {v, x") (A. 3) 



--S JJSJx^'" 

is not. If the interfacial thickness is small compared to the radii of curvature, the difference between (f> and (pa is 
small and vanishes for planar interface considered under cartesian coordinates. However, it is $ but not which is 
a physical amount, and thus is the surface density. 

2. On the integration of a gradient function in curvilinear coordinates. 

a. Consider a function (j) being the divergence of a vector function: (j) = V-q(r). Then 
(V-q)'=^(r;x") = V-q(r) - {V-qff(r)} e(a;'' - xi) - {V-q'^(r)} e(a;i - x") 

(A.4) 

= V-(q'=^)(r;a;'') -fqf(r)-Ve(a;'' -xi) -t-q^(r)-Ve(2:i - x") 
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where 

q'=^(r; a;'*) = q(r) - q5(r) 6(2;" - Xi) - q^(r) e(xi - x") (A.5) 
and q^ and q^ are defined similarly to cj)^ and (j>^. Furthermore its excess 

^(^''1-11) = tAi / dxi/ii/i2/»3(V-qr(r;a:«) (A.6) 

Using the standard formula for the divergence of a vectorial function in curvilinear coordinates 

^ / d d d \ 

V-q^ , , , i T^{h2h3qi) + T^ihih3q2) + T^{hih2q3)) (A.7) 

hi /i2 /13 \dxi dx2 0x3 J 



one can show that 



[ dxihih2h3V-iq^^){v-x'') = 

+ dxi {—{hih3qr) + j^ihih2qr)] 

g,s J^g.s ydx2 ^8x3 ) 



= h2h3q1 



(A., 



Using the standard formula for the gradient of a scalar function in curvilinear coordinates 

V6I-— — i +— — i (A 9) 

hi dxi /i2 dx2 /13 dxs 

one can show that for Heaviside step function 



dxi/ii/i2/i3q''(r)-Ve(a;i-a;") = q''(2;^r||)•il = gl(a;^r,l) (A.IO) 

Substituting Eq. (|A.4p into Eq. (|A.6|) and using Eq. (|A.8D and Eq. (|A.10|) we obtain 

V^(a;^r||) = gi(a;^r||)-gi(a;^r||) + V^|(a:^r|,) (A.ll) 

where 

^1 - hk C (^^('^^'^3.2-) + ^('^I'^.C)) (A.12) 

and we took into account that according to Eq. (jlll.ip q'j^ix^''') = q'j^{x^''') = 0. 

b. Consider a special form of a vector q for which q|| = J||(/) where V|| -Jn = 0. Here V|| is a parallel component 
of three-dimensional nabla-operator so that 

^11-^11 =h;h^3{^2^'''''"''^^i;^^'^"''^ 

(A.13) 

1 de . 1 86 



Then 



h2 0x2 h3 0x3 



-g^ihih3q2) + -g^ihih2q3) = hMs (J||-(V||0) + (V||-J||)0) = /ii/i2/i3 J|| • (V||(/)) (A.14) 
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